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Abstract 

In this study, the hnear, third-order nonhnear and total refractive index (RI) 
changes of mult i- layered quantum dot (MLQD) and multi-layered quantum 
anti-dot (MLQAD) with a hydrogenic impurity are calculated. Our numer- 
ical results indicate that the RI changes are considerably sensitive to the 
geometrical parameters of systems and the incident optical intensity. It is 
observed that by changing the incident photon energy, the RI curves cor- 
responding to MLQAD and MLQD are considerably different in shapes and 
behaviors. Generally, our results show that MLQAD models have very higher 
RIs than MLQD models. 

Keywords: Quantum dots and anti-dots. Optical intensity. Refractive 
index. Confining potential 



1. Introduction 

Spherical-shaped semiconductor nanostructures such as quantum dots 
(QDs) and quantum anti-dots (QADs), have gathered a great deal of at- 
tention in the last decades [l[ [i], y, 0, la, 0]. These confined structures, have 
atomic-like discrete energy levels (subbands) and particular optical proper- 
ties. Recently, investigations of the physical properties of some new nanos- 
tructures such as the multi-layered quantum dots (MLQDs) and anti-dots 
(MLQADs) have attracted the attention of researchers 0,0, y. These new- 
found spherical nano-systems can be created by synthesizing different semi- 
conductor materials. The investigation of properties of such systems when 
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they are doped with shallow donor impurities is one of the most studied 
problems during the last decade , [m, [12I, [1^ . Impurity existence in such 
systems plays the essential role in their electrical and optical properties. The 
nonlinear optical properties (absorbtion coefficient and refractive index) as- 
sociated with optical absorption of mentioned nanostructures are known to 
be much stronger than bulk material, due to the quantum confinement ef- 



fect 14, 15, 16 



17, 18, 19 



These particular optical properties have the 
highly advantage and potential to construct the optoelectronic devices, such 
as photo-detector, quantum dot laser and high speed electro-optical modu- 



lators [20|, |2l|, |22|, |23| . Due to relevance of QDs and QADs to several tech- 
nological applications, their linear and nonlinear absorption coefficient and 
refractive index changes have been investigated both theoretically and ex- 



perimentally by many authors [2J, |25|, l26|, l27|, |28|, l29|, |30|, iMl, |32|, |33 



The optical properties of QDs was considered by Efros [2^ for the ffist 
time. He studied the light absorption coefficient in a spherical QD with in- 
finitely height walls. In ref 26|, the authors have calculated not only the lin- 
ear and nonlinear absorption coefficients but also the refractive index changes 
in a three-dimensional Cartesian coordinate quantum box with finite confin- 
ing potential barrier height. The numerical simulations in 3J| show the RI 
changes in a parabolic cylinder QD and indicate that RI changes are strongly 
affected by the dot size, the optical intensity and the polarization of electro- 
magnetic field. The linear, nonlinear and total refractive index changes and 
absorption coefficients for transitions in a spherical QD with parabolic poten- 
tial have recently been studied in ref 32|. Their results, expressed in several 
allowed transitions Is — Ip, Ip — Id and Id — If, show that the transition 
between orbitals with big / (orbital quantum number) values move to lower 
incident photon energy region in the presence of parabolic potential term. 



The outline of this paper is as follows. The theory and formulation for 
both MLQAD and MLQD are briefly presented in Section 2. The linear, 
nonlinear and total RIs of both models are plotted for various conditions 
as the function of incident photon energy in Section 3. In this Section the 
comparison between the behaviour of both models are presented. A brief 
summary is presented in the last section. 
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H Ga^ ^Al^As with potential V| 



H Ga^ ^Al^.As with potential 



Figure 1: A schematic view of the MLQD and MLQAD structures 



2. Theory and formulation 

2.1. Hamiltonian and Schrddinger equation 

By uniting of GaAs, Gai-xAl^As and Gai„j,AlyAs materials, one can 
make the multi-layered nanostructures. We introduce two types of these 
adjacent material connections next to each other as shown in Fig. [H More 
precisely, a MLQD (MLQAD) consists of a spherical core made of GaAs 
(Gai_j^AlyAs) surrounded by a spherical shell of Gai_a;Al^.As (Gai_2,.Al^.As), 
encompassed in the bulk of Gai_j^Alj^As (GaAs). In other words, the MLQAD 
is made whenever the core material and the bulk material in MLQD change 
places. We denote a for the core radius and b for the total dot (core plus 
shell) radius, therefore, shell thickness is 6 — a. 

Let us consider an electron confined in the MLQD (MLQAD). Within 
the framework of the effective mass approximation, the Hamiltonian of the 
system, with hydrogenic impurity in center, is given by 

fc2 2 

H = ^V^-^ + Wr) (1) 

where /i(r) and e are the electronic effective mass and dielectric constant in 
the semiconductor medium. In the different regions, m(r) is represented as 
follows 

r < a 

m{r) = /i2 a < r < b (2) 

r > b 
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V{r) is the confining potential (CP) that corresponding to MLQAD and 
MLQD has a following form 



{Vq r < a 

a<r<b (3) 
V2 r>b 

and 

{V2 r < a 

Vi a<r<b (4) 

Vo r>b 

The barrier height Vi arises due to a mismatch between the electronic affini- 
ties of the adjacent media in MLQAD and MLQD, so Vi (V2) corresponds to 
binding the Gai_a;Alj.As (Gai_yAlj,As) to GaAs (Gai_a;Ala,As). The value of 
Vq is always zero (i.e., the GaAs material has no CP) but writing Vq enable 
us to write next equations in the compact form as you will see. 

For many QD heterostructures, such as GaAs/Gai_a;A^As, the image 
charge effects and polarization can be notable in the MLQD and MLQAD 
if there is a large dielectric discontinuity between the core dot and the sur- 
rounding medium. However, this is not the case for the GaAs/Gai_a;Ala;As 



multi-layered systems [38[, therefore these effects can be ignored safely in our 
calculations. Furthermore, for the sake of generality, the difference between 
the electronic effective mass in the dot (antidot) core, the shell and the bulk 
materials have been ignored (i.e., /xi = /i2 = /^s = )• The energies are mea- 
sured in meV, the effective Rydberg, Ry = /ie^/2^^(47re)^, and distances 
are expressed in = 4:7ieh? /{fie"^). For instance, in the particular case of 
GaAs-based semiconductors, /i = 0.067me, and e = 13.18eo- Thus, for a 
GaAs host the effective Rydberg and the effective Bohr radius numerically 
are Ry = 5.2meV and oq = 10.4nm respectively. 

The Schrodinger equation is Hijj{r, 9, ip) = Etp{r, 9, ip), by use of separa- 
tion variable method and this fact that the hydrogenic potential is spherically 
symmetric, the eigenfunction in the spherical coordinate can be rewritten as 



ip{r,9,ip) = R{r)Y{9,(p). Y{9,(p) is given by the spherical harmonics [36 
and it is independent of radial components thus we only focus on the radial 
part, R{r), when we investigate the optical properties of system. The radial 
part {R{r)) is affected not only by the core and total radius but also by the 
CPs. 
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2.2. The radial Schrddinger equation for MLQAD and MLQD 

The radial part of Schrodinger equation, in the spherical coordinate, for 
both MLQAD and MLQD model takes the form 

fif, +-^- ^^)R^\r) + {E-V, + ^)R^\r) = (5) 
IjjL dr"^ r dr r^ Arcer 

where i = 2, 1 and (z = 0, 1 and 2) are corresponding to r < a, a < r < b 



and r > b for MLQAD (MLQD) respectively. Recently in we have 
analytically solved this equation for both MLQAD and MLQD and we have 
shown that the solutions of radial Schrodinger equation can be expressed 



in terms of Whittaker and hypergeometric functions (One can refer to [11 
for reading with full detail) which could be very complicated, or almost 
impossible to be applied to the MLQAD or MLQD to satisfy the boundary 
conditions. Therefore, we apply the finite difference method to calculate the 
direct numerical solutions for Eq. (5). 

2.3. Linear, nonlinear and total refractive index 

The linear and the third-order nonlinear optical RI changes for the inter- 
subband transitions can be calculated by the density matrix approach and 
the perturbation expansion method. For this purpose the system under study 
can be excited by an external electromagnetic field of frequency u, such as 
the following 

E{t) = Ee'""' + Ee-''^' (6) 

If E is the perpendicular electromagnetic field along the z axis, the Hamilto- 
nian of system becomes Hq + ezE{t) where where Hq is the Hamiltonian of 
system without the electromagnetic field E{t). The linear and the third-order 
nonlinear optical RI of a spherical nano-system, within a two-level system, 
in a special case, fro m g round (ipwo) to first allowed excited state (V'210) can 



be expressed as [26|, |29 



(7) 
(8) 



rir 2nfeo {E21 - fkoy + {hV 21)^ 

An(3)(a;) _ fiIca^\M2i\^ {E21 - hu) 

Ur ~ ^0 [(^21 - + (/ir2l)2]2 

Then the total RI is 

= I (9j 

^yij-y- T'lj'Y* t^tTp 
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In above equation we have used 



E21 — E2 — El, 

M21 = I < '02io|e^|'0ioo > I, (10) 

where £'21 denotes the difference of the energy between two lowest electronic 
states, M21 is an element of electric dipole moment matrix, /x, cr^ and rir 
are the permeability, carrier density and refractive index of the system re- 
spectively, huj is the incident photon energy, r2i = I/T21 where T21 is the 
phenomenological relaxation rate, caused by the electron-phonon, electron- 
electron and other collision processes. / is the optical intensity of incident 
wave, c and £0 is the speed of light in the free space and permeability of 
vacuum. 



3. Results and discussion 

In the following, we calculate the linear, third-order nonlinear and to- 
tal refractive index changes in both nano-systems with different shell thick- 
nesses and CPs. The unchanged parameters used in our calculations are: 
r2i = 0.2ps, (7^ = 3.0 X 1022^-^ Ur = 3.2, £r = 13.18. 

According to Eq. (7) and (8), the linear and the nonlinear RI always 
have an opposite sign and their values depend on the value of E21. More 
exactly, in tkij — E21 both terms are zero and change their sign as tkij < E21 
{huj > E21) the linear term is positive (negative) and the nonlinear term is 
negative (positive). Furthermore, the smaller £"21 leads to the larger linear 
and nonlinear values. In other words, the smaller £21 causes the larger value 
of RI terms and also for smaller E21, the RI terms reach zero value in smaller 
incident photon energy. 

Under condition {/ = 200MW/m'^, a = 2ao, 6 - a = 0.5ao, Vi = 2Ry 
and V2 = 5Ry}, our numerical calculations result that the values of E21 for 
nano-systems are fj^^Q^^ = Q,37rneV and E^^^*^^ = lQ.98meV, thus we 
have 

^MLQAD ^ ^MLQD ^^^^ 

According to discussion in above discussions, one can conclude that the value 
of linear and nonlinear RI for MLQAD should be very larger than correspond- 
ing value in MLQD and also linear and nonlinear RI for MLQAD should reach 
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Figure 2: The linear, nonlinear and total RIs as a function of incident photon energy in 
the same condition {/ = 2QQMW/m'^, b - a = 0.5ao, Vi = 2Ry and V2 = 5Ry} for (a) 
MLQD (b) MLQAD 

zero value in very smaller incident photon energy than MLQD as Fig.[2]shows 
well. These figures indicate that in MLQAD (MLQD) the major and minor 
terms of total RI are nonlinear (linear) and linear (nonlinear) terms respec- 
tively. 

From relation (7)- (9), by increasing the optical intensity the major term 
(nonlinear term) of total RI of MLQAD increases although the minor term 
remains unchange, so we predict that the total RI for MLQAD should be 
enhanced. But in the case of MLQD, by increasing the optical intensity the 
major term (linear term) of total RI do not change although the minor term 
increases. Since the opposite sign of major and minor terms, we predict that 
by increasing the optical intensity the total RI should be decreased. Fig. [3] 
verifies our predictions. 

In Figs, m the total RIs of MLQAD are plotted when the shell thickness 
is fixed as 6 — a = Ioq whiles a core radii have a four different values as 
a = 2, 2.5, 3, S.Soq. In Figs. S] (a) and (b) the values of CPs are {Vi = 
2Ry, V2 = 5Ry} and {Vi = 2Ry, V2 = ooRy} respectively. It is observed 
that by increasing the core dot radius, the position of total RI peak has no 
significant shift but the height of total RI peak enhanced. Furthermore, by 
changing the value of V2, no considerable differences are observed. 

The same conditions are employed for MLQD in Figs. El In this case, 
by increasing the core dot radius, the maximum value of RI becomes bigger 
and also the peak position shifts toward the lower photon energies. This 
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Figure 3: The linear, nonlinear and total RIs as a function of incident photon energy with 
the fixed parameters {a — 2ao, b — a = 0.5ao, Vi = 2Ry and V2 = 5Ry} and three different 
intensities for (a) MLQD that from up (b) MLQAD that from bottom the intensity value 
is / = 200MVK/m2, AOOMW/m'^ and 600MW/m'^ respectively. 




Figure 4: The total RI of MLQAD as a function of incident photon energy with the fixed 
shell thickness b — a — Igq and four different core radii for (a) Vi — 2Ry and V2 — 5Ry 
(b) Vi = 2Ry and V2 = ooRy 
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Figure 5: The total RI of MLQD as a function of incident photon energy with the fixed 
sheU thickness b — a = lag and four different core radii for (a) Vi = 2Ry and V2 = 5Ry 
(b) Vi = 2Ry and V2 = ooRy 



behaviour arises from this fact that by increasing the core dot radius, the 
energy difference between ground and first excited states becomes smaller. 
But in the case of MLQAD (Fig. H]) the energy difference is approximately 
negligible, thus the peak position remained unchanged. Simultaneously com- 
parison of Fig. |5] (a) and (b) shows that RI value not affected by changing 
in the second potential barrier (V2). 

The total RI curves of both nano-systems for the fixed core radius a = 
2.0ao and four different shell thicknesses b — a = 1, 0.5, 0.25, 0.02ao, as the 
function of photon energy are presented in Fig. [6] and Fig. [71 respectively. 
The difference between (a) and (b) parts in these figures is the value of 
V2 that is 5Ry and infinity respectively. In Fig. |6] and Fig. [71 one can see 
that by increasing the shell thickness, the RI value corresponding to both 
nano-systems becomes greater. It is observed that by increasing the shell 
thickness, the peak position of total RI curves for MLQAD approximately 
remains unchange, however in MLQD the peak and valley of total RI shift to 
smaller photon energies. The influence of second CP (V2) on total RI values 
is more prominent in the case of different shell thickness rather than constant 
shell thickness. In Fig. [6] and Fig. [71 comparing the curves for different V2 
within certain shell thickness for both nano-system show that the changes 
corresponding to smaller shell thickness is more clear rather than others. 
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Figure 6: The total RI of MLQAD as a function of incident photon energy with a — 2ao 
and four different sheh thicknesses for (a) Vi — 2Ry and V2 = 5Ry (b) Vi — 2Ry and 
V2 — ooRy 
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Figure 7: The total RI of MLQD as a function of incident photon energy with a — 2ao 
and four different shell thicknesses for (a) Vi ~ 2Ry and V2 = 5Ry (b) Vi — 2Ry and 
V2 = ooRy 
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4. Conclusion 

We have investigated the hnear, the third-order nonhnear and the to- 
tal optical RIs of multi-layered spherical nano-systems with a hydrogenic 
impurity in the center. As our results indicate, the optical RIs changes of 
MLQAD and MLQD are very different. It is found that in the same condi- 
tions MLQAD has very higher Rl changes than MLQD. It is observed that 
for the fixed shell thickness, by increasing the amount of core radius, the 
total RI peak heights become larger and shift toward lower incident photon 
energies for MLQD but the peak heights corresponding to MLQAD become 
larger with no significant changes in peaks positions. Also, in this case the, 
changes in V2 value lead to no considerable changes in total RI curves. Fur- 
thermore, for a fixed core radius, by increasing the shell thickness, the peak 
of total RI curves corresponding to MLQAD, moves to larger values without 
considerable changes in the incident photon energy but the peak heights re- 
lated to MLQD become larger and shift to smaller energy regions. Finally, it 
is found that in the latter case, by decreasing the shell thicknesses, the total 
RI curves are considerably affected by V2 value. 

References 

[1] D. Hciss, M. Kroutvar, J. J. Finley, G. Abstreiter, Solid State Commun. 
135 (2005) 591. 

[2] M. Grundmann, D. Bimbcrg, N.N. Ledentsov, Quantum Dot Het- 
erostructures, Wiley, Chichester, 1998. 

[3] Y. Masumoto, T. Takagahara, Semiconductor Quantum Dots, Springer, 
Berlin, 2002. 

[4] Y. R Varshni, Phys. Lett. A 252 (1999)248. 

[5] E. Sadeghi, Physica E 41 (2009)1319. 

[6] B. Vaseghi, G. Rezaei, V. Azizi, F. Taghizadeh, Physica E 43 (2011) 
1080. 

[7] C. Y. Hsieh, D.S. Chuu, Phys. J Condens. Matter 12 (2000) 5641. 
[8] Y. Naimi, A. R. Jafari, J Comput Electron. 11 (2012) 414. 



11 



[9] M. J. Karimi, G. Rezaei, J. Appl. Phys. Ill (2012) 064313. 
[10] H. Ta§, M. §ahin J. Appl. Phys. Ill, (2012) 083702. 
[11] A. R. Jafari, Y. Naimi, J Comput Electron. 

[12] S. Davatolhagh, A. R. Jafari, M.R.K. Vahdani, Superlattices and Mi- 
crostruchure. 51 (2012) 62. 

[13] V. A. Holovatsky, O.M. Makhanets, O.M. Voitsekhivska, Physica E 41 
(2009) 1522. 

[14] R. Atanasov, F. Bassani, V. M. Agranovich, Phys. Rev. B 50 (1994) 
7809. 

[15] L. Tsang, D. Ahn, S.L. Chuang, App. Phys. Lett. A 362 (2007) 37. 

[16] M. J. Karimi, G. Rezaei, Physica B 406 (2011) 4423. 

[17] W. Yang, X. Song, R. Li, Z. Xu: Phys. Rev. A 78 (2008) 023836. 

[18] W. Xie, S. Liang, Physica B 406 (2011) 4657. 

[19] U. Bockelman, G. Bastard, Phys. Rev. B 45 (1992) 1688. 

[20] M. Troccoh, A. Belyanin, F. Capasso, E. Kubukcu, D. J. Sinco, A.Y. 
Cho, Nature 433 (2005) 845. 

[21] X. Jiang, S. S. Li, M. Z. Tidrow, Physica E 5 (1999) 27. 

[22] C. Liu, B. Xu, Phys. Lett. A 372 (2008) 888. 

[23] C. S. Lent, P.D. Tougaw, W. Porod, G.H. Bernstein, Nanotechnology 4 
(1993) 49. 

[24] A. Al. Efros, A. Efros, Sov. Phys. Semicond. 16 (1982) 772. 

[25] G. Rezaei, B. Vaseghi , F. Taghizadeh , M. R. K. Vahdani, M.J. Karimi, 
Superlattice and Microstructures 48 (2012) 450. 

[26] Unlu, S., Karabulut, I., and Safak, H.: Physica E 33 (2006) 319. 

[27] W. Xie, Optics Communications 283 (2010) 3703. 



12 



[28] X. C. Li, AM. Wang, Z. L. Wang, Y. Yang, Superlattice and Microstruc- 
tures 51 (2012) 580. 

[29] C. Zhang, Z. Wang, M. Gu, Y. Liu, K. Guo, Physica B 405 (2010) 4366. 

[30] S. Liang, W. Xie, H. A. Sarkisyan, A. V. Meliksetyan, H. 
Shen, Superlattice and Microstructures 51 (2012) 868. 

[31] I. Karabulut, S.J. Baskoutas, App. Phys. 103 (2008) 073512. 

[32] B. Cakir, Y. Yakar, A. Ozmen, Journal of Luminescence 132 (2012) 
2659. 

[33] A. Ozmen, Y. Yakar, B. Cakir, U. Atav, Optic Communications 282 
(2009) 3999. 

[34] M. R. K. Vahdani, G. Rezaei, Phys. Lett. A 374 (2010) 637. 
[35] R. Kordad, Indian. J. Phis 86 (2012) 513. 

[36] N. Zettih, Quantum Mechanic, Concept and Application, Wiley, New- 
York, 2001 pp. 287-289. 

[37] A. Abramowitz, I. Stegun, Handbook of Mathematical Function with 
Formulas, Graphs and Mathematical Tables, US GPO, Washington 
(1994) pp. 505-509. 

[38] S. Adachi, J. Appl. Phys. 58 (1985)R1. 



13 



